Accepted for publication in the Astronomical Journal - 2009 
Mar 19 



Spitzer Observations of Comet 67P/Churyumov-Gerasimenko at 

5.5-4.3 AU From the Sun 

Michael S. Kelley/'^ Diane H. Wooden,^ Ceciha Tubiana^, Hermann Boehnhardt,^ Charles 

E. Woodward,^ David E. Marker,^ 

ABSTRACT 

We report Spitzer Space Telescope observations of comet 67P/Churyumov- 
Gerasimenko at 5.5 and 4.3 AU from the Sun, post-aphehon. Comet 67P is the 
primary target of the European Space Agency's Rosetta mission. The Rosetta 
spacecraft will rendezvous with the nucleus at heliocentric distances similar to our 
observations. Rotationally resolved observations at 8 and 24 //m (at a heliocentric 
distance, Vh, of 4.8 AU) that sample the size and color-temperature of the nucleus 
are combined with aphelion i?-band light curves observed at the Very Large 
Telescope (VLT) and yield a mean effective radius of 2.04 ±0.11 km, and an 
i?-band geometric albedo of 0.054 ± 0.006. The amplitudes of the i?-band and 
mid-infrared light curves agree, which suggests that the variability is dominated 
by the shape of the nucleus. We also detect the dust trail of the comet at 4.8 
and 5.5 AU, constrain the grain sizes to be < 6 mm, and estimate the impact 
hazard to Rosetta. We find no evidence for recently ejected dust in our images. 
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If the activity of 67P is consistent from orbit to orbit, then we may expect the 
Rosetta spacecraft will return images of an inactive or weakly active nucleus as 
it rendezvous with the comet at = 4 AU in 2014. 

Subject headings: Comets: individual (67P/Churyumov-Gerasimenko) — In- 
frared: solar system — meteoroids 



1. INTRODUCTION 



Ecliptic comet 67P/Churyumov-Gerasimenko (67P) is the primary target of the Rosetta 
mission, which will rendezvous with the nucleus in 2014 when the comet has passed aphelion 
and is 4 AU from the Sun (aphelion distance 5.71 AU). Rosetta will orbit with the comet 
through perihelion in 2015 (perihelion distance 1.25 AU). The Rosetta spacecraft is designed 
to characterize the comet nucleus (morphology, composition) and coma (development of 
activity, dust-gas interaction, interaction with the solar wind), with both orbiting and landing 



spacecraft (Glassmeier et al. 2007). Understanding the comet's gas and dust behavior near 



aphelion is important for mission planning. 

In previous perihelion passages, the comet was not characterized over the orbital arc 



at which the Rosetta spacecraft will rendezvous with the nucleus (Agarwal et al. 2007b). 



Dynamical models by FuUe et al. (2004) of the dust coma and tail of 67P observed near 



perihelion in 2002-2003 have indicated that the comet is significantly active and able to 
eject millimeter and centimeter sized grains at 3.6 AU, post-aphelion (i.e., pre-perihelion) . 



However, this conclusion is disputed by Agarwal et al. (2007b), who point out that the 



dust tail analyzed by FuUe et al. (2004) contains overlaping contributions from all dust 



ejected before the comet reached 1.5 AU. Direct imaging of the comet with ground-based 



observations have shown that the comet is active at 2.0 AU, post-aphelion (Agarwal et 



al. 2007b). Recent ground-based optical and Spitzer Space Telescope mid-infrared (mid- 



IR) observations have extended that limit to 2.9-3.0 AU (Kadota et al. 2008; Wooden et 



al. 2008). Some ecliptic comets are active at more than 4 AU from the Sun, e.g., comets 



65P/Gunn, 74P/Smirnova-Chernykh, 152P/Helin-Lawrence, and 103P/Hartley (Lowry & 
Fitzsimmons [2001 : [Lowry et al. 1999; Snodgrass et al. 2008). Indeed, even comet 67P was 
brighter than expected for a bare nucleus and dust trail in Spitzer spectra of the comet at a 
heliocentric distance (r^) of 4.98 AU, pre-aphelion (Kelley et al. 2006), and in ground-based 
CCD observations at 4.9 AU, pre-aphelion (Mueller 1992 ; Lamy et al. 2006 ) — equivalent post- 



aphelion observations have only recently been presented (Tubiana et al. 2008). The present 



study of comet 67P is motivated by the lack of aphelion observations and the importance of 
the characterization of the nucleus and comet activity of 67P for the Rosetta mission. 
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In this paper, we present and discuss Spitzer Space Telescope (Gehrz et al. 12007 Werner 



et al.|2004 ) observations of comet 67P at heliocentric distances between 5.5 and 4.3 AU, post- 



aphehon — the latter distance is comparable to the Vh at which Rosetta will rendezvous with 
the comet nucleus. In §2, we present imaging observations of the comet obtained with both 
the Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004) at 24 fim and the 
Infrared Array Camera (IRAC; Fazio et al. 2004) at 8 /xm, and 14-35 /im spectra of the 
nucleus with the Infrared Spectrograph (IRS; Houck et al. 2004). In 24 /im MIPS images, 
we detect the comet's dust trail (Sykes & Walker 1992) and a tail composed of large grains. 
We constrain the largest grain size ejected by 67P to be < 6 mm, and then estimate the 
large grain number density and its impact hazard to Rosetta in Our observations are 

consistent with an inactive nucleus between 5.5 and 4.3 AU, post-aphelion. We present 
upper-limits to emission from dust at 5.5 and 4.8 AU in ^3.2 In ^ 3.3[ we derive the infrared 
(IR) light curve at 4.8 AU by sampling the emission from the nucleus at 11 epochs with 
consecutive MIPS and IRAC observations. The Spitzer derived light curves of comet 67P 
constrain the shape, color-temperature, and effective radius of the nucleus. We examine 
those properties and their behavior with rotational phase to show in ^3.4 that the visible 
light curve is likely dominated by nucleus shape rather than albedo. However, there are 
some discrepancies between the visible and IR light curves, therefore some albedo variation 
cannot be ruled out. Our results are summarized in ^ 



2. OBSERVATIONS AND REDUCTION 
2.1. Imaging 

Spitzer observed comet 67P with the MIPS 24 /im instrument in Photometry/Super 
Resolution mode (pixel scale 2.55 arcsec pixel"^). This mode provides 14 dithered images 
covering an approximate field-of-view of 7' x 8' (each image is approximately 5' x 5'). The 
comet nucleus is centered in the dither pattern, and the telescope tracks the comet with 
its computed non-sidereal rates. Each image has an exposure time of 10 s. The dither 
pattern was executed 9 times for each of 2 epochs at = 5.5 AU (2006 August 16 and 
2006 September 01), 3 times for each of 11 epochs at 4.8 AU (2007 May 18), and once 
for each of 3 epochs at 4.4 AU (2007 September 19). The comet is readily identified as a 
moving point source, centered in the field-of-view, and bisected by the comet's dust emission 
(Fig. [T]), except at 4.4 AU where dust is not observed due to crowding by point sources 
and diffuse emission arising from the interstellar medium. For the 5.5 and 4.8 AU epochs, 
we median combined each observation's images and divided each image with the result to 

— the 4.4 AU images are too crowed for 



remove image artifacts (Spitzer Science Center|[2008 
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this step. We stacked the images together in the rest frame of the comet with the MOPEX 



software (Makovoz & Khan 2005) to produce 2 mosaics at 5.5 AU, 11 at 4.8 AU, and 3 at 

4.4 AU. The 11 MIPS 24 fim images at 4.8 AU spanned a duration of 13.17 hr (start to 
finish) enabhng samphng of the hght curve over approximately one rotational period of the 
nucleus {Prot ~ 13 hr, see ^3.3). Table [T] presents the observation parameters for the MIPS 
data. 

[Table [l|, Figure [l) 

We also observed comet 67P with the IRAC 8 fim camera in mapping mode (pixel scale 
1.22 arcsec pixel"^). The 11 observations spanned a duration of 13.09 hr (start to finish) 
and, similar to the MIPS observations, measured the light curve just over one rotational 
period of the nucleus at rh = 4.8 AU (2007 May 17). The observations were designed with a 
14-point cyclic dither pattern, and 30 s exposure times. The images of each observation are 
mosaicked in the rest frame of the comet with the MOPEX software. The nucleus is faint at 
8 /im and not readily identifiable. We median combined all 11 images and found the comet 
to be near the center of the frame (Fig. [T]). The comet was not detected with the IRAC 

4.5 fim camera, as expected for a i? = 2 km nucleus at this distance. Table [T] presents the 
observation parameters for the IRAC data. 

Because the IRAC and MIPS instruments are not simultaneously operated, Spitzer 
observed the two light curves as close in time as was operationally feasible. The IRAC 
observations were executed first, one rotational period was skipped for instrument calibration 
and changeover, then the MIPS observations were obtained. We assumed a rotational period 
of 12.8-12.9 hr to constrain each light curve using approximately equal time steps between 
each observation. This allows us to phase the two wavelengths (8 and 24 /im) together 



without interpolation. Tubiana et al. (2008) recently constrained the rotational period to be 
12.7047 ±0.0011 hr at aphelion. Using this new period, the phase offsets between individual 
IRAC and MIPS observations are minimal and range from 0.005 to 0.022 rotational periods, 
with a median offset of 0.014. Our light curve sampling of 67P differs from the approach of 



Lamy et al. (2008) who observed the comet at 24 fim with Spitzer 16 times — each observation 



was separated by 0.5, 1.0, or 1.5 hr, spanning a total of ~ 12.5 hr. 



2.2. Spectra 

We obtained two spectra of comet 67P with the IRS long wavelength (14-35 yum), low- 
resolution (A/AA ~ 100) module. One spectrum was taken while 67P was 4.9 AU from the 
Sun, and the other at 4.3 AU from the Sun. The parameters of the IRS observations are 
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presented in Table [TJ The default IRS Stare observing template was used, which placed the 
comet at two positions separated by 55" in the 11" wide by 168" long slit. Spitzer positioned 
the slit at the comet's predicted ephemeris coordinates, and followed the comet with its com- 
puted non-sidereal rates. The observations were repeated 9 times for a total of 18 spectra. 
We differenced the two slit positions to remove the background, then replaced bad pixels 
with the average value of the nearest neighbors or ignored them altogether. We extracted 
spectra from the images using the Spitzer IRS Custom Extraction (SPICE) software pack- 
ag^ The profile of the comet is equivalent to a point source, therefore we utilized the default 
point source aperture provided by SPICE, which is 7.2 x A/27.0 fim pixels wide (pixel scale 
5.1" pixel"^). The individual spectra are averaged together to produce two final spectra, one 
at each of 4.9 AU and 4.3 AU. The spectra, smoothed with a 5-point statistically weighted 
average, are presented in Fig. [2j 

[Figure |2] 



3. RESULTS AND DISCUSSION 
3.1. The Dust Environment 

3.1.1. An Empirical Model 

Understanding the near-nucleus dust environment of 67P when Rosetta will rendezvous 
with the comet is important for mission planning. Our 24 /xm images contain emission from 
large dust grains. With photometry of this dust, we can estimate the number density of 
large grains near the nucleus, and assess their impact hazard to spacecraft. Furthermore, if 
we can adequately describe the morphology of the dust thermal emission it can be removed 
from the images prior to the measurement of the comet nucleus. To remove the dust, we 
fit surface brightness profiles to linear cuts (3 x 51 pixels) measured perpendicular to the 
projected velocity vector of the comet, then subtract our best-fit model from the images. The 
assumed profile is the sum of a Gaussian function (the dust) and a linear term (the image 
background). The results of our empirical model fits to the = 4.8 and 5.5 AU images are 
presented in Figure [3] and Table [2] The photometry has been color-corrected (x 1.057) to 
the MIPS 24 yum central wavelength of 23.7 /im using a blackbody spectral template with 
Td = 300 r^°'^ K, where is the typical dust effective temperature of trail grains (W. T. 



Reach et al., 2009, in preparation; Sykes & Walker 1992). Two trends are evident in the fits 



^SPICE is available at http : //ssc . spitzer . caltech . edu/ 
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to the dust surface brightness: 1) the peak surface brightness hnearly varies with distance 
from the nucleus, and 2) within the uncertainties (± a few degrees in position angle) the dust 
emission is aligned with the orbit, as would be expected for a dust trail. To generate model 
images of the dust, we approximate the morphology with the best-fit parameters described 
in Table [2j The model dust images are subtracted from the 5.5 and 4.8 AU MIPS images. 

[Table g Figure |3] 



3.1.2. A Dynamical Model 

Previous images of 67P at aphelion show a neck-line tail structure. A neck-line structure 
is a thin, hnear feature caused by the projected alignment of dust grains that were emitted 



at a true anomaly 180° prior to the observation in question (Southworth 1963 Kimura & 
Liu||1975" ). The enhancement is strongest when the observer is located near the orbital plane 
of the comet, which is frequently the case for ecliptic comets such as 67P. Small dust grains 
cannot remain near the nucleus for many years; therefore, if a neck-line is observed in our 
Spitzer observations, it must be comprised of large dust grains that weakly interact with solar 
radiation pressure. Observations of neck-lines are beneficial because they allow us to measure 
the size distribution of dust grains emitted by the comet. This becomes apparent when we 
consider that larger dust grains are ejected from the nucleus with lower ejection velocities and 
are accelerated less effectively by solar radiation pressure than smaller grains. Therefore, the 
largest grains remain the closest to the nucleus, and the grain size decreases with increasing 
distance along the neck-line. In principle, the variation of the surface brightness distribution 
of a neck-line with distance from the nucleus can be transformed into a grain size distribution. 
In the Spitzer MIPS images, the mere existence of a neck-line can help us constrain the largest 
grains emitted by 67P. The age of a neck-line observed on 2006 August 16 would be 3.9 years, 
corresponding to dust emitted at a post-perihelion distance of = 1.26 AU. Similarly, if a 
neck- line were observed on 2007 May 18, its age would be 4.6 years, i.e., released from the 
nucleus at rh = 1.30 AU, post-perihelion. 

In order to determine if a neck-line tail structure exists in our images, we compute 



simulated images of the comet using the dynamical dust model of Kelley et al. (2008). The 
primary forces acting on a dust grain are the force from solar radiation (F^), the force of 
solar gravity (Fg), and gravitational perturbations from the planets. Dynamical models of 
comets typically describe dust grains with the parameter (3 = Fr/Fg. Both forces vary with 
r^^, therefore (3 is independent of heliocentric distance. The [3 parameter does, however. 
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depend on grain size: 



/5 



0.57Q 



pr 



i<g pa 

where Qpr is the radiation pressure efficiency (~ 1 for large grains), p is the grain density in 
units of g cm~^, and a is grain radius in units of pm. (Burns et ah 1979). When converting 



from /?- values to grain radii, we follow Kelley et al. (2008) and adopt a grain density of 
1 g cm~^ (see Kelley et al. (2008) for discussion related to this choice). To describe the dust 



ejection dynamics and history of comet 67P, we start with the best parameters of Kelley et 



al. (2008). In summary, Kelley et al. (2008) fit Spitzer MIPS images of the comet taken on 



2004 February 23 (4.5 AU, pre-aphelion) and found three distinct dynamical components: 
1) a dust trail with an age greater than one orbit, 2) a neck-line ejected near the 2002 
perihelion passage, and 3) a broad dust tail. The surface brightness distribution of the dust 
in the 4.5 AU, pre-aphelion, image was best modeled by: 

• a dust production proportional to cos Zq, where Zq is the sun-zenith angle of a surface 
element ejecting dust grains. 



a size distribution proportional to a 



-3.5 



an ejection velocity equal to 0.5 j3^'^ ' cos 2;© km s ^, and 
dust grains with —5 < log/? < — 1 (6 pm < a < 60 mm). 



Kelley et al. (2008) assumed the comet's dust production is proportional to dust production 
rates inferred from optical photometry Kidger (2003), Qd oc r^^'^. We note that the above 



parameters are similar to those of Ishiguro (2008), derived by comparison of dynamical mod- 



els to three epochs of ground-based optical observations. The greatest differences between 



the Kelley et al. (2008) and Ishiguro (2008) models is the choice of the largest grain sizes 

-5.8 



sizes (60 mm versus 5 mm) 



and dust production rates {Qd oc versus 



We generated two simulations of 67P using the dynamical model and best parameters 



of Kelley et al. (2008), with the added stipulation that no grains are ejected outside of 
r/i = 4.5 AU. This distance corresponds to the comet-Sun distance during the February 



2004 Spitzer observations, where little to no coma was observed (Kelley et al. 2008 Lamy 



et al. 2008). The simulations reproduce the observing conditions of our 5.5 AU and 4.8 AU 
images. Each simulation tracked ^ 2 x 10^ grains and recorded their final positions, which 
are projected onto the celestial sphere for an observer located at Spitzer. The oldest grains 
tracked were ejected near aphelion in March 1986, and during the simulation, the comet 
passed perihelion three times. The resulting images are convolved with a model 24 pm point 
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spread function generated with version 2.0 of the Tiny Tim/ Spitzer softwar^ In Fig. |4] and 
[5} we present the simulated 24 /im images and surface brightness profiles, generated with 
a variety of /5-value cutoffs: (a) /5 > 1 x 10-^ (b) /5 > 5 x 10-^ (c) /? > 1 x 10-^ (d) 
/? > 2 X 10""^, (e) /5 > 5 X 10~^. Simulations (c) and (d) best reproduce the observed profile 
at both epochs. We estimate that 67P ejects grains with /3-values larger than 1 x 10~^, 
corresponding to grains smaller than ^ 6 mm in radius. This is essentially the same size 



upper-limit derived by Ishiguro (2008). 



Formally, a neck-line tail structure is present in our Spitzer observations. The neck-line 
and dust trail overlap in both the 5.5 AU and 4.8 AU images, and they cannot be resolved into 



two components as has been observed by Kelley et al. (2008), Ishiguro (2008), and Agarwal 



et al. (2007a). The distinction between neck-line and trail in our observations rests merely 



on their dynamical parameters, since their grain populations have approximately the same 
sizes (~ 1 mm in radius) at these heliocentric distances. Taking into account our /5-value 
lower-limit of 10~^, the neck-line is located at a position angle ~ 1° north (at 5.5 AU) or 
south (at 4.8 AU) of the trail. In our images, we find that 15-25% of the observed flux within 
6.5" of 67P's projected orbit is attributable to the neck-line tail, released at ~ 1.3 AU. 



3.1.3. Large Dust Grain Number Density and Impact Hazard to Rosetta 



Kelley et al. (2008) estimate the trail surface brightness near the nucleus to be ~ 0.1 



MJy sr ^ in 24 /im MIPS images obtained at 4.5 AU pre-aphelion. Their estimate depends 



on a dynamical model to separate trail emission from tail emission. Kelley et al. (2008) then 
derive a number density of ~ 10~^^ m~^ for millimeter sized grains. From our empirical 
model (Table [2]), the peak dust surface brightness at the nucleus is 0.076 ± 0.004 MJy sr~^ 
and the FWHM is (19±2) x 10^ km at 4.8 AU, post-aphelion. These values are improvements 



over the Kelley et al. (2008) analysis because the dust width and surface brightness may be 



measured closer to the nucleus. Moreover, we can estimate the extent of the dust along 
the Spitzer line-of-sight using the above dynamical simulations. Figure |6] shows the 4.8 AU 
simulation as viewed by an observer located above the ecliptic plane. The apparent width 
of the dust along Spitzer's line-of-sight is ~ 60, 000 km. 

[Figure |6] 

Assuming the dust is entirely comprised of millimeter sized grains with a dust tempera- 
ture of Td = 300 rr°-5 K and contained within a volume of 38, 000 x 60, 000 x 40, 000 km, we 



■^Tiny Tim/ Spitzer is available from the Spitzer Science Center: 



http://ssc.spitzer.caltech.edu/ 
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derive a grain number density of (1.33 ±0.06) x 10~ m~ . Our assumed volume corresponds 
to the observed width of the trail (2cr, Table [2]), the thickness of the simulated dust along the 
line-of-sight in the (3 > 10~'^ model, and the MIPS spatial resolution. The number density 



is about a factor of 10 smaller than the estimate of Kelley et al. (2008). We caution that 



the uncertainty is only based on the observed flux, and we are assuming that our chosen 
simulation parameters accurately represent the behavior of the comet. 

Inspection of Fig. [6] reveals that a spacecraft, such as Rosetta, approaching from the 
interior of the comet's orbit will not encounter the bulk of the dust. To estimate an upper- 
limit to Rosetta^s impact hazard based on our simulations and observations, we assume the 
dust grains are uniformly distributed in space, and the Rosetta^s path length through that 
dust is 30,000 km (half of the dust's thickness in comet's orbital plane). The total impact 
probability to millimeter-sized grains is < 0.003% m~^. Taking the spacecraft cross-section 



to be 4 m^ and the total solar panel cross-section to be 64 m^ (Glassmeier et al. 2007), 



Rosetta has < 0.3% chance of encountering a millimeter-sized grain during approach to the 
nucleus. During comet fly- by missions, spacecraft typically survive collisions with a few 0.1- 
1 mg (?a 500 /im) dust grains at high relative speeds, ~ 1 km s~^ ( A'Hearn et al.||2005 Green 
et al. 2004 McDonnell et al. 1987 1993). However, previous fly-by missions visited their 



respective comets during vigorous coma activity, which enhances the density of large grains. 
Therefore, our encounter probability upper-limit is not unusual and seems reasonable. 



3.2. Comet Activity 

The Spitzer MIPS observations are suitable for a search for faint coma emission. We 
restricted our search to the MIPS images because they are more sensitive to emission from 
dust at 4-6 AU from the Sun. First, the 11 MIPS observations at 4.8 AU are averaged 
together to create a deep image, then the empirical dust model is subtracted from the 
images. Subtracting the dust model should not affect our search for a dust coma as comae 
are peaked near the nucleus, yet the observed dust varies linearly across the position of the 
nucleus (Fig.jsj). The nucleus of 67P has a full- width at half-maximum (FWHM) of 2.2 pixels, 
equivalent to the FWHM of point sources in the surrounding field-of-view when the images 
are mosaicked in the celestial reference frame. A comparison of a model PSF generated 
with Tiny Tim/ Spitzer (T=200 K) to the observed radial profile of 67P (Fig. [T]) strongly 
suggests that the nucleus is devoid of an extensive coma at this heliocentric distance. The 
azimuthally averaged coma flux at 6-8 pixels (15-20") from the nucleus, i.e., outside of the 
first diffraction ring of a point source, is 0.0025 ± 0.0030 MJy sr~^, yielding a Icr upper-limit 
of 0.0055 MJy sr~^. Similarly, we also measured a stellar profile and a coma upper-limit of 
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0.0042 MJy sr"Mn the images at 5.5 AU. 
[Figure [T] 

We can convert the surface brightness upper-hmits, /(p), to an integrated coma bright- 
ness upper-hmit if we assume a coma shape. We assume a steady-state coma surface bright- 
ness distribution of p~^, where p is the angular distance from the nucleus. This assumption 



is typical for photometry of distant comets (e.g., Snodgrass et al. 2008). The more steeply 



shaped profile of a radiation pressure dominated coma, I^, oc p~ -^■^ ( |JewittfcMeech|1987D , has 
a less restrictive upper-limit by a factor of 2. The integrated coma brightness upper-limits 
at 24 pm are computed from Fc < 27rp^/(p), and are equal to 0.26 mJy at = 4.8 AU, 
and 0.20 mJy at 5.5 AU. These upper-limits are 9% of the point source fluxes presented in 
^ 3.3 Removing a 9% coma contribution from our photometry decreases our best-fit nucleus 
radii by 19%. In ^3^, we assume the coma contribution is zero, and find that the resulting 
nucleus radii agree with previous estimates. 



3.3. The Mid-infrared Light Curve and the Nucleus Size, Shape, and 

Temperature 

3.3.1. Photometry 

The MIPS and IRAC light curves allow us to constrain the surface temperature dis- 
tribution of the nucleus, which is required in order to accurately model its size and shape. 
First, we describe our techniques for measuring the comet flux in our light curve observations 
at 4.8 AU; next, we describe our techniques for the snapshot observations at 4.4 and 5.5 AU; 
and finally, we describe our nucleus thermal model and the results of the model fits to the 
data. 

For the images taken at 4.8 AU, we subtract the empirically derived model dust from 
each mosaicked MIPS observation. Next, we subtract the celestial background. As the 
telescope tracks the comet, background objects approach and recede from the vicinity of the 
nucleus. To remove the celestial objects, we identify image pairs where the comet has moved 
at least 6 pixels. Aligning and differencing these image pairs removes the fixed background 
objects. Finally with aperture photometry, we derive the flux of the comet utilizing for IRAC 
a 5 pixel source aperture, and a 10-25 pixel background annulus. The aperture correction for 
the IRAC photometry is x 1.082, and the color-correction is x 0.869. For MIPS the aperture 
and color-corrections are, respectively, x 1.180 — 1.349, and x 1.047 for 4-6 pixel source 
apertures and a 15-30 pixel background aperture. The color-corrections were generated and 
iteratively refined using a best-fit near-Earth asteroid thermal model (NEATM; Harris||1998 ) 
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spectrum (see below), as prescribed by the Spitzer Science Center (2006, 2008) 



At 4.4 and 5.5 AU, a slightly different approach is required. First, our empirical dust 
model is subtracted from the 5.5 AU images. At 4.4 AU, the images were not sufficiently 
sensitive to detect a dust trail, therefore no trail was subtracted. The nucleus of 67P falls 
very close to a few background stars in all images. The comet was not tracked sufficiently 
long enough to allow celestial sources to move away from the nucleus, therefore the shift 
and align technique we used to remove point source contamination in the 4.8 AU images 



is not possible. Instead, we subtract these stars with a PSF fitting program (Diolaiti et 



al. 2000). The program does not preserve the critically-sampled Spitzer PSF adequately 
enough for precise photometry (it is designed for over-sampled images), but does sufficiently 
remove background stellar sources fainter than the nucleus for our purposes. We repeat the 
same photometry procedure as executed on the 4.8 AU data (color-corrections are x 1.046 
at 5.5 AU and x 1.050 at 4.4 AU). 

Table [T] presents the color- and aperture-corrected photometry. The light curve photom- 
etry has also been corrected to remove the effects of the changing observation geometry as 
the distance between the Earth and 67P increased during the 1.5 days that elapsed between 
the first IRAC and the last MIPS observation (< 1.5% correction). The absolute calibra- 
tion errors for MIPS and IRAC are 3% and 4%, respectively, but are not included in the 
reduced photometry in order to enable direct comparison of the point-to-point variation of 
the data. Absolute calibration errors are included in the thermal model fit results discussed 
below. The mean 24 /im and mean 8 /im fluxes at 4.8 AU are F24 = 2.83 ± 0.15 mJy and 
Fg = 0.109 ±0.015 mJy. 



3.3.2. Nucleus Thermal Model 



To derive 67P's nucleus size and temperature, we adopt the NEATM (Harris 1998) 



noting that this model assumes the nucleus is spherical and in instantaneous equilibrium 
with insolation. Modeling the nuclei of comets with the NEATM has been successful in the 
past ( Fernandez et al.|2006 Kelley et al.||2006 Fernandez et al.pOOO ), including a successful 
reproduction of the global properties of the nucleus of comet 9P/Tempel substantiated by 



resolved temperature maps created from Deep Impact observations (Groussin et al. 2007). 



Despite the simplifications of the model, the NEATM can still provide good estimates for 
the effective radius and albedo of the comet (e.g., see Groussin et al. 2007). 



In the NEATM, the temperature distribution of the surface peaks at the sub-solar point 
and decreases with the cosine of the Sun-zenith angle. The hemisphere beyond the solar 
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terminator has a temperature of K. The sub-solar temperature is computed via 



T 



1 



r]ea 



(K), 



(2) 



where A is the Bond albedo, S = 1365r^^ W m~^ is the solar flux at the distance of the 
comet {vfi measured in AU), t] is the IR-beaming parameter, e = 0.9 is the infrared emissivity, 
and a is the Stefan-Boltzmann constant. Commonly, the geometric albedo, p, is reported 
in the literature, rather than the Bond albedo {p = A/q, where q is the phase integral, 
see Manner et al. 1981 Bowell et al. 1989). The range of geometric albedos derived for 
comet nuclei is 2-6% ( Lamy et al.||2004 ), which has a minor impact on surface temperature. 
7] remains as the primary unknown quantity. Formally, rj parameterizes how the thermal 
emission is preferentially emitted toward the Sun due to surface roughness, but because 
the NEATM assumes instantaneous equilibrium with insolation, t] must also account for 
the thermal inertia of real world surfaces. In practice, t] is constrained by the shape of the 
nucleus thermal emission. We derive t] by iteratively comparing the observed spectral energy 
distribution to the spectrum of a NEATM modeled hemisphere, observed at the same phase 
angle. The primary source of uncertainty in rj is our IRAC 8 /im photometry. 

Since the thermal emission from comet nuclei is weakly dependent on the surface albedo. 



we turn to optical measurements to constrain this value. We follow the methods of Rus- 



sell (1916) and Jewitt (1991), to employ the relationship between radius, absolute visual 



magnitude, and albedo: 



p 



1.496 X 10« 
R 



1 2 



^q0.2{Mq-H) 



(km). 



(3) 



where R is the effective radius in km, p is the geometric albedo, H is the absolute optical 
magnitude of the nucleus, and the Mq is the absolute magnitude of the Sun. Using a 
reduced-x^ fitting routine, we simultaneously vary p, rj, and R in Eq. [3]and the NEATM to 
constrain these parameters for comet 67P. 

The non-simultaneous IRAC and MIPS light curves must be phased together using a 
measurement of the nucleus rotation period. The cadence and time span of our Spitzer 
observations do not permit unique determination of the rotation period of the nucleus. 



However, other investigators have more appropriate observations. Lamy et al. (2006) ob 



served 67P with the Hubble Space Telescope in March 2003, and constrained the period to 
P = 12.41 ±0.41 hr. Lowry et al. (2006) observed the comet over three nights with the NTT 



in May 2005, and measured P = 12.72 ±0.05 hr. Tubiana et al. (2008) have phased together 



VLT light curves measured in May 2006 and August 2006 (near aphelion), and determined 
the period to be 12. 7047± 0.0011 hr, commensurate with the previously determined periods. 
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although with greater accuracy. Because of the high accuracy, and that it was obtained clos- 
est in time to our observations, we adopt the aphelion period of 12.7047 hr for our analysis 
of the Spitzer data. 

The variation of the nucleus' color-temperature with rotational phase can be assessed 
by combining the IRAC and MIPS light curves at 4.8 AU, and fitting a Planck function to 
each pair of data points. The derived color-temperatures are within la of the mean value, 
187 ± 4 K, at all longitudes. Thus, we do not detect any variation in temperature with 
rotation period. 



As discussed in [:3.4 the aphelion light curves from Tubiana et al. (2008) agree in 



amplitude with our Spitzer observations. Using our mean fluxes and their mean i?-band 
absolute magnitude {Hr = 15.35 ± 0.04, Mq^r = -27.14), we derive R = 2.04 ± 0.11 km, 
rj = 0.68 ± 0.06, and pr = 0.054 ± 0.006, via minimization. We note that our fit has no 
degrees of freedom because we are fitting three data points with a three parameter model. 
The mean effective radius agrees with other recent measurements, which range from 1.9 
to 2.1 km once differences in albedo, absolute solar magnitude, and adopted phase angle 



behaviors are taken into account (Lamy et al. 2007 Tubiana et al. 2008 Tubiana 2008). 



The axial ratio, derived from the MIPS light curve minimum (2.382 mJy) and maximum 



(3.37 mJy), is a/b > 1.41 ± 0.07, which also agrees with the analyses of Lamy et al. (2006) 



and Tubiana et al. (2008). The low value strongly suggests the surface has a low thermal 



inertia, similar to the surfaces of other comets ( Groussin et al.||2007 ). The best- fit radii from 
all our MIPS observations, derived assuming r] = 0.68 ± 0.06, are presented in Table [ij 



3.4. Comparison to Other Light Curves 



Lamy et al. (2008 ) created a 24 /im light curve derived from observations of 67P obtained 



when the comet was at 4.5 AU from the Sun, pre-aphelion. This latter light curve is compared 



to our data in Fig. 9 We scaled the Lamy et al. (2008) data by: 



NEATM(rfe,i,A,,i,a,,i,?/,p^,i?) 
NEATM(r;,,2,A,,2,a.,2,r/,Pi?,i?) 



0.746, 



(4) 



where is the comet- Spitzer distance, as is the Sun-comet- Spitzer angle, subscript 1 desig- 
nates the observing conditions for our MIPS light curve (Table [T]), subscript 2 designates the 
observing conditions for the Lamy et al. (2008) light curve (r/j = 4.48 AU, = 4.04 AU, 
as = 12.1°), rj = 0.68, and pr = 0.054. The scale factor is independent of R. The absolute 
phasing of the two light curves is unknown — over 2000 rotational periods have elapsed be- 



tween the two observations. We have shifted the Lamy et al. (2008) light curve to minimize 
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the point-to-point variation between the two data sets. Only near 0.1 phases do the two 
hght curves significantly disagree, but since both light curves are sparsely sampled at these 



phases the discrepancies cannot be resolved. Lamy et al. (2008) combine their mid-IR light 



curve with a shape model of the nucleus, derived from visible data, and they estimate the 
i?-band albedo lies in the range 0.039-0.043. This value agrees with our best-fit value at the 
2cr level, and we attribute the discrepancy to a difference in the adopted values of the comet's 
absolute optical magnitude. Indeed, if we use the absolute magnitude Hji = 15.46 derived 
by Lamy et al. (2008), our derivation of the comet's albedo decreases to pn = 0.049 ± 0.006, 



which coincides within la of the Lamy et al. (2008) estimates. We prefer to adopt the op- 



tical magnitude from Tubiana et al. (2008) since it was obtained when the comet was near 



aphelion, similar to our Spitzer observations. 
[Figures |9) [To) 



Observations of 67P's visible light curve (Tubiana et al. 2008 Tubiana 2008) indicate 



that the overall variability profile is the same throughout the aphelion arc of the cometary 
orbit. Hence, we take the well-sampled -R-band light curve of 67P from the May and Au- 
gust 2006 Very Large Telescope (VLT) data as a reference for comparison with the Spitzer 



measurements (Fig. 10). Visual inspection of the light curves reveals that the maximum-to- 



minimum flux ratios (1.41 and 1.45) are similar. This result suggests that the variability at 
visible through mid-infrared wavelengths is dominated by the effective radius of the nucleus 
as it rotates, rather than variations in the albedo. 

Since it is not possible to do an exact phasing of the infrared and visible light curves, we 
applied arbitrary shifts to the visible light curve in order to force the maxima into agreement 



(Fig. 10, top). The shifted light curves have significant differences at 0.1-0.6 phases, and are 
systematically different at 0.4 to 0.6 phases. We also attempted to minimize the point-to- 



point variations between the mid-IR and i?-band light curves (Fig. 10, bottom). The overall 



agreement has improved in this case, but significant differences still exist. The disagreements 
may indicate deviations in scattering and/or thermal behaviors of the surface, e.g., variations 
in albedo, surface roughness, or thermal conductivity. Or, there may be systematic errors in 
the photometry, which could explain a few rogue data points. For example, moving sources 
in the background are not accounted for with our background removal techniques, but may 
be faint and not apparent in the images. Such background sources have the greatest effect 



on the faintest parts of the light curve. It is interesting to note that Lamy et al. (2008) also 



had some difficulty finding a perfect agreement between their mid-IR light curve and a shape 



model derived using a light curve inversion technique on optical photometry (Lamy et al. 



2007). A further analysis of these discrepancies would require detailed shape, thermophysical. 



and scattered light models, better constrained phase shifts, and/or simultaneously obtained 



mid-IR and visible light curves. If we consider the overall agreements in peak-to-peak fluxes 
and the approximate agreement in light curve shapes, then we may place confidence in our 
nucleus albedo, color-temperature, and effective radius computed in §3.3.2 



3.5. Nucleus Spectra 

In principle, our two spectra of the nucleus can provide an independent measurement 
of the nucleus effective-temperature. Both spectra (Fig. [2]) were taken during the post- 
aphelion arc (r/j = 4.9 and 4.3 AU) along with our MIPS and IRAC observations. Using 
the NEATM, we fit two thermal models to each unbinned spectrum. In one model the 
IR-beaming parameter is a free parameter, and in the other model it is fixed at t] = 0.68 
( ^3.3.2[ ). For both models the albedo is fixed at pr = 0.054. Our model results are presented 
in Fig. |2| and the radii derived from the constant rj model fits are listed in Table [1} The 



absolute calibration error for the IRS instrument is approximately 10% (Spitzer Science 



Center 2008), which is included in our radius uncertainties. 



The IRS spectrum at Vh = 4.9 AU does not significantly constrain the IR-beaming 
parameter [rj = 1.3 it 0.4). Furthermore, when the IR-beaming parameter is allowed to vary, 
the best-fit radius increases to 3.0 ± 0.5 km to compensate for the high best-fit //-value. 
Such a high ?7-value is unphysical given the constraints from our imaging data. A more 
meaningful IR-beaming parameter and effective radius is derived from the 4.3 AU spectrum: 
R = 1.99 ± 0.21 km, and rj = 0.80 ± 0.14. This 77-value agrees with the value of 0.68 ± 0.06 
derived from the MIPS-IRAC light curve. Reduced values range from 0.6-0.8 for all 
model fits (?7-fixed and r/- variable, 4.9 AU and 4.3 AU). 



4. SUMMARY 



We observed comet 67P/Churyumov-Gerasimenko between 5.5 and 4.3 AU from the 
Sun, post-aphelion, including a full light curve of the nucleus at 8 and 24 fim at 4.8 AU. 
We employ the NEATM to derive a mean effective radius of R = 2.04 ± 0.11 km, and an 
average IR-beaming parameter of r/ = 0.68 ± 0.06. There is no evidence for variations of 
the surface temperature with rotational phase. The amplitude of the light curve constrains 
the primary-to-secondary axis ratio to a/b > 1.41 it 0.07. The mid-IR light curve is in 



good agreement with the MIPS 24 /xm light curve presented by Lamy et al. (2008), when 



both datasets are phased with a 12.7047 hr rotation period ( Tubiana et al.|[2008 ), and when 



the Lamy et al. (2008) light curve is phase shifted to minimize point-to-point variations. 
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Together, the composite mid-IR hght curve roughly agrees with two aphehon i?-band hght 



curves measured by Tubiana et al. (2008) in May and August 2006. Significant discrepancies 



exist between the two data sets, which may indicate surface variations of albedo, surface 
roughness, or thermal conductivity. However, the maximum-to-minimum fiux ratios agree, 
suggesting that the nucleus variability is dominated by the shape of the nucleus as it rotates. 
Using the absolute magnitude 
i?-band geometric albedo to be 



15.35 ± 0.04 (Tubiana et al. 2008), we estimate the 
0.054 ±0.007. 



The Spitzer spectrophotometric observations of the central point source are wholly 
commensurate with the thermal emission from an inactive nucleus. Therefore, we conclude 
that vigorous coma activity did not occur before = 4.3 AU in 2007. As the Rosetta 
spacecraft approaches the comet during the same in-bound path in 2014, we may expect to 
see images of an inactive, or weakly active nucleus. 

We also observed the dust trail and large dust grain tail (i.e., neck-line tail structure) 
surrounding the nucleus in 24 /xm images taken at 4.8 and 5.5 AU from the Sun. Our 
other images were not sensitive enough to detect the dust. Starting with the dust ejection 



parameters of Kelley et al. (2008) and Ishiguro (2008), and additionally imposing a maximum 



grain size of 6 mm {(3 > 10 ^), simulations of the aphelion dust environment are consistent 
with our MIPS 24 yum observations. Our maximum grain size agrees with the maximum 



grain size of 5 mm derived by Ishiguro (2008). The grain number density near the nucleus 
is (1.33 ± 0.03) X 10~^^ m~^, assuming the dust is composed of millimeter-sized grains. The 
density corresponds to a Rosetta dust impact probability upper-limit of 0.3% during the 
spacecraft's approach to the nucleus in 2014. 
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Fig. 1. — Spitzer images of comet 67P/Churyumov-Gerasimcnko. The nucleus is located at 
the center of all images and marked with vertical bars for clarity. The image scale is provided 
in (a). North, east, the projected Sun angle (0), and the projected comet velocity vector 
(v) are indicated in the lower- right of each image: (a) a co-added mosaic of all 11 IRAC 
8 iim observations from 2007 May 17 (4.8 AU); (b) a co-added mosaic of all 11 MIPS 24 /im 
images from 2007 May 18 (4.8 AU); (c) the MIPS 24 //m image taken on 2006 September 01 
(5.5 AU); and (d) a MIPS 24 //m image from 2007 September 19 (4.4 AU). The dust trail 
crosses two of the MIPS images (b and c) along the velocity vector. 



- 24 - 



10 
8 



6 r 

4 - 




8 - 



-I — I — I- 



Variable rj 
Fixed rj 



2007 Apr 1 9 (4.9 AU) 



-| — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — h 

2007 Sep 30 (4.3 AU) 




15 



20 25 30 

Wavelength {jurni) 



35 



Fig. 2. — Spitzer IRS 14-38 jim. spectra (circles) of comet 67P taken at rh = 4.9 and 4.3 AU. 
The spectral resolution has been reduced with a 5-point statistically weighted average. Two 
best-fit thermal models are shown for each spectrum, one with a variable IR-beaming pa- 
rameter (77), and one with the IR-beaming parameter fixed at 0.68 ± 0.06. At 4.9 AU: 
R = 3.0± 0.5, r] = 1.3± 0.4, xl = 0-62 (dashed-line), and R = 2.11 ± 0.17, r] = 0.68 ± 0.06, 
Xl = 0.78 {soM-line). At 4.3 AU, R = 1.99 ±0.19, r] = 0.80 ±0.14, xl = 0.62 (dashed-line) , 
and R = 1.82 ± 0.14, r] = 0.68 ± 0.06, xl = 0-78 {solid-line). 
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Offset along trail (") 



Fig. 3. — Peak trail surface brightness and best-fit linear trends derived from our empirical 
trail model fits (Table ^ for three MIPS 24 fxm epochs. Positive x-axis values lie along the 
projected velocity vector (i.e., these positions are ahead of the comet nucleus). The data 
points in the grayed area were excluded from the profile fits to prevent the nucleus from 
affecting the results. 
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Fig. 4. — Top: Simulated 24 /xm images (a-e) and the 2006 Aug 16 MIPS image (f). Each 
simulated image was generated with a different large grain cutoff: (a) /? > 1 x 10"^, (b) 
/3 > 5 X 10-^ (c) /3 > 1 X IQ-^ (d) /3 > 2 X IQ-^ (e) /3 > 5 x IQ-^ All images are scaled 
to a value of 1.0 at a distance of 2' behind the nucleus. North is up, east is to the left, the 
projected Sun angle (©) and comet velocity vector (v) are indicated with arrows, and the 
image scale is provided in (a). Bottom: Surface brightness profiles generated from each of 
the simulated images (a-c) and the observation (f). We fit a linear profile to the observed 
dust in (f), and plot it as a dotted-line in each subplot. Vertical lines mark the location of 
the nucleus. 
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Fig. 5. — Same as Fig. |4], except the 2007 May 18 epoch simulations and MIPS image are 
shown. 
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(a) S3 (b) 




Fig. 6. — Simulated 24 /xm images of 67P's 2007 May 18 dust environment as viewed by 
an observer located above the ecliptic plane: (a) all dust with (3 > 10~^; and (b) neck-hne 
tail dust with (3 > 10~^. A solid-line marks the projected orbit of the comet, vertical lines 
mark the location of the nucleus, the projected Sun angle (©), comet velocity vector (v), 
and Spitzer direction (S) arc indicated with arrows, and the image scale is provided in (a). 
The surface brightness of image (b) has been scaled by a factor of 2 relative to image (a) 
to improve contrast, and both images have been convolved with a 4600 km wide (FWHM) 
Gaussian kernel. The neck-line tail grains and other dust ejected near perihelion have orbits 
that are distinctly different from the trail grains, which are located closer to the orbit of the 
comet. 
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Fig. 7. — The 24 fim radial surface brightness profile (triangles) of comet 67P/Churyumov- 
Gerasimenko at 4.8 AU (2007 May 18) after emission from the dust trail has been removed. 
A 200 K model point spread function is plotted with dots, and the la background level is 
marked with a horizontal dashed-line. There is no evidence for a dust coma in this or any 
other image. 
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Fig. 8. — Spitzer MIPS 24 //m (triangles) and IRAC 8 //m (squares) light curves, phased 
with a rotational period of 12.7047 hr. The IRAC light curve has been scaled by a factor of 
10. 
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Fig. 9. — Spitzer MIPS 24 /im (triangles) light curve (this work), and the pre-aphehon 



(rh = 4.5 AU) MIPS 24 /xm light curve from Lamy et al. (2008) (diamonds). We have 
phased the light curves with a 12.7047 hr rotational period. The pre-aphelion light curve 
has been shifted in phase space and scaled by 0.78 to qualitatively match our post-aphehon 
data. 
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Fig. 10. — Same as Fig. lol but now the May and August 2006 i?-band light curves (crosses] 



from Tubiana et al. (2008) have been added. The i?-band light curve has been shifted to 



qualitatively match the Spitzer data. Top: The i?-band data has been shifted to match the 
mid-IR light curve maximum. Bottom: The i?-band data has been shifted by A(f)r = —0.13 
phases so that the light curves better agree at 0.3-1.0 phases. 



